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Abstract. Let p be a prime number. We present an algorithm that computes 
p-integral bases in number fields. The algorithm is obtained as a by-product 
of a p-adic factorization method based on Newton polygons of higher order. 
The running-time of the algorithm appears to be very good; it computes the 
2-integral basis of a number field of degree 1152 in a few seconds. 

1. Introduction 

The theory of higher order Newton polygons developed in [Mon99j and revised 
in [HN| (HN standing for "higher Newton" ) has revealed itself as a powerful tool 
for the analysis of the decomposition of a prime p in a number field. Higher order 
Newton polygons are p-adic objects, and their computation involves no extension 
of the ground field, but only extensions of the residue field; thus, they constitute an 
excellent computational tool. Based on this technique, in [GMN08] we presented an 
algorithm, due to the second author, that computes the prime ideal decomposition 
of p in a number field, and the p- valuation of the discriminant. In this paper we 
show that a slight modification of this algorithm provides a p-integral basis too, as 
a by-product. 

Let f{x) e Z[x] be a monic irreducible polynomial of degree n. Let 9 & Q 
be a root of f{x), K = Q(0) the number field determined by /(x), and "Lk the 
ring of integers of K. Along the flow of Montes' algorithm [GMNOSl Sects. 2, 3], 
some special polynomials (t){x) are constructed such that the </>-adic developments 
of f{x) provide the combinatorial data that are necessary to build the higher order 
Newton polygons of f{x). Just by keeping trace of the quotients that are obtained 
along the computation of these ^-adic developments, we can construct n integral 
elements of the form q{9)/p^ , where q{x) is a polynomial that is a product of some 
of these quotients, and the exponent v of the denominator is expressed in terms 
of combinatorial data attached to the polygons (Proposition 13. 6p . We conjecture 
that these elements are always a p-integral basis. The conjecture is proven in two 
extreme (and opposite) cases: when the algorithm uses only polygons of order one, 
regardless of the number of p-adic irreducible factors of f{x) |EN09| . and when 
f{x) is irreducible over Zp, regardless of the orders of the Newton polygons that 
are considered by the algorithm (Theorem 14. 4[) . Thus, the remaining task is to 
prove the conjecture in all in-the-middle (or mixed) cases. 

In spite of the apparent uncompleteness of the paper, the algorithm is quite 
useful in practice because it does compute a p-integral basis. In fact, the p- valuation 
ind(/) v{1jk'- ^[6*]) of the index of f{x) is one of the outputs of the algorithm; 
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therefore, it is trivial to check numericaUy that the output family of n integral 
elements is a p-integral basis indeed. More specifically, if we denote by O the Z- 
module generated by this output family, the algorithm computes the p-valuation 
of {Zfc : O) and prints this value together with a message of "maximal order" or 
"non-maximal order". Thus, the algorithm either provides a p-integral basis or it 
warns the user that this were not the case. On the other hand, in an exhaustive 
numerical exploration that we carried out, we always obtained the output message 
of "maximal order" . 

This method is very efficient and it makes possible to compute p-integral bases 
in number fields of high degree. In fact, the running-time of the algorithm appears 
to be very good. Even in some bad case, chosen to test the limit of its capabilities, 
it computes the p-integral basis of a number field of degree 1920 and p-index 958560 
in a few seconds, in a personal computer. However, the p-integral bases produced 
by the algorithm are not triangular; if we add at the end of the algorithm a trian- 
gulation process, the running-time may increase in a significant way, because the 
matrix of the coefficients of the polynomials q{x) tends to be far from sparse. One 
can try to avoid this by considering only generators of I^K.p '■— "^k ®z ^(p) as a 
Z(p) [0]-module and by using a Groebner Basis algorithm over to find a p-stem 
of the integral basis (a set of generators of I^K.-p as a Z(p) [0]-module, with minimal 
degrees). This procedure leads to a much faster triangulation process, but at the 
cost of an increase of the memory requirements. 

The outline of the paper is as follows. In section 2 we review Montes' algorithm. 
In section 3 we study the partial quotients qix) obtained along the computation of 
0-adic developments of /(x), and we determine the highest exponent v such that 
q{d)/p'^ are integral fProposition l3.6|) . In section 4 we construct our candidates to 
be a p-integral basis and we describe how to adapt Montes' algorithm to compute 
them as a by-product. In Theorem 14.41 we prove that these integral elements are a 
p-integral basis when fix) is irreducible over Zp[a;]. In section 5 we present some 
numerical explorations with times of execution of the algorithm in comparison with 
the standard procedures to get a p-integral basis of PARI, MAGMA and SAGE. 

Conventions and general notations. We fix a number field K ~ Q{d), 
generated by a monic irreducible polynomial f{x) G Z[a;] of degree n, such that 
f{9) — 0. We denote by Zk the ring of integers of K. 

We fix also a prime number p E Z. We denote by v. Qp — > Q, the p-adic 
valuation normalized by w(p) = 1; we extend w to a discrete valuation of Qp{x) by 
letting it act as 



on polynomials. We shall introduce p-adic valuations Vr of higher order of the same 
field Qp{x), and we convene that vi = v. 
We denote by 



the p-adic value of the index of the polynomial f{x). Recall the well-known rela- 
tionship, i!(disc(/)) = 2ind(/) + v{disc{K)), between ind(/), the discriminant of 
/(x) and the discriminant of K. 




ind(/) ■.= v{{Zk: Z[0])), 
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If F is a finite field and <fi{y), ipiv) G ^[y], we write (p ^ ip to indicate that the 
two polynomials coincide up to multiplication by a nonzero constant in F. 

We shall freely use the notations and results of |HN| . A short review can be 
found in ^GMNOSl Sect.2]. 

2. Review of Montes' algorithm 

The flow of Montes' algorithm can be represented by a tree. The starting node 
is labelled by the input polynomial /(x), and the first branches that sprout from 
this node are labelled by the different irreducible factors of f{x) modulo the input 
prime p. For each one of these branches, we lift the corresponding irreducible factor 
modulo p to some monic polynomial (t>i{x) £ Z[a;]; then, we compute the (/)i-adic 
development of f{x): 

f{x) = aa{x) + ai{x)(l)i{x) H h as{x)cl)i{xy , 

with ai{x) £ "^[x], degai{x) < deg0i(x). This development leads to a Newton 
polygon of the first order, iVi(/), and we compute its principal part, N^{f), that 
gathers only the sides of negative slope. The length of this polygon N^{f) (that is, 
the length of its projection to the horizontal axis) is the highest exponent ni such 
that divides f{x) modulo p. 

Caution: there is a different Newton polygon of the first order for each branch, 
but we use the notation (/) whenever it is implicitly understood that we are 
working on a concrete branch. 

To each side of (/), with slope A, we attach a residual polynomial, R\{ f)(i/) £ 
Fi[y], with coefficients in the finite field Fi := Z[x]/{p, (l)i{x)). The branch we are 
analyzing ramifies into as many branches (of first order) as pairs (Xjiply)), where 
A is one of the slopes of N^{f) and ?A(y) is a monic irreducible factor of R\{f){y) 
in ¥i[y]. 

Let us fix now one of these first order branches, and let {4>i{x); Ai, ^i{y)) be the 
data that characterize the branch. Once these data are fixed, for any polynomial 
P[x) £ Z[x] we shall denote simply by Ri{P){y) the residual polynomial attached 
to the side of slope Ai of the Newton polygon Ni{P) constructed from the 4>i- 
development of P{x). We construct now a polynomial 4>2(x) such that: 

(1) 4>2{x) is congruent to a power of 0i(x) modulo p, 

(2) Ni{(j)2) is one-sided, with slope Ai, 

(3) i?i(02) ~^i(2/) inFi[2/]. 

The data {(f>i{x); Xi,ipi{y)) determine also a p-adic valuation V2 of the field 
Qp{x). We can use this valuation of second order to construct Newton polygons of 
second order of any polynomial P{x) G with respect to its 02-development; 

if P{x) = '^Q<:iai{x)(l>2{xy , then N2{P) is the lower convex envelope of the set 
of points {i,V2{ai{x)(j)2{xy)) of the Euclidean plane. The algorithm computes the 
principal part N2 {f)', to each side of this polygon, with slope (say) A, we attach a 
residual polynomial of second order, R\{f){y) G F2[?/], with coefficients in the finite 
field F2 := Vi[y]/{ij)i{y)). The length of N2{f) is equal to the highest exponent n2 
such that "01(2/)"^ divides Ri{f){y). The branch we are analyzing ramifies into as 
many branches of second order as pairs {X,ip{y)), where A is one of the slopes of 
^2 if) i'iy) is a monic irreducible factor of R\{f ){y) in F2[y]. 
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Going on with this procedure, a branch of order r of the algorithm is labelled by 

a type of order r: 

t = (01 (x); Xi,(j)2{x)\ . . .]Xr-l,cj)r{x)\ Xr,i>r{y)), 

where (/>i(a;), . . . , 4>r{x) £ Z,[a;] are monic polynomials, Ai, . . . , Ar G Q~ are negative 
rational numbers, and ipr{y) G ^r[y] is a monic irreducible polynomial. With these 
data one is able to construct ap-adic valuation Vr+i and a special monic polynomial 
4>r+i{x), called the representative of the type, and the pair Wr+i, (t>r+iix) allows us 
to build (r + l)-th order Newton polygons. 

For the precise definition of a type see [HNl Sect. 2.1]. For our purposes, it is 
sufficient to recall that for all 1 < i < r, these objects satisfy the following pro- 
perties: 

(1) The Newton polygon of i-th order of (pi^i{x) with respect to t, Ni{ipij^i), 
is one-sided, with slope A^. 

(2) The residual polynomial of «-th order of with respect to t satisfies: 
Ri{4>i+i){y) ^ i^i{y) in Fi[y], for certain monic irreducible factor 4'i{y) of 

in F, [£/]._' 

(3) The length of N^_^_^{f) is the highest exponent n^+i such that 
divides R^{f){y) in F,[y]. 

Every type carries implicitly a certain amount of extra data, whose notation we 
fix now. For all 1 < i < r: 

• hi, Ci are a pair of positive coprime integers such that Ai — —hi/ei. 

• fi := degi/jiiy), and /o := toi := deg0i(x). 

• mi+i := deg4)i+i{x). Note that m^+i = eJ^mi = /o(ei/i) • • ■ (eifi). 
Also, for all 1 < i < r -t- 1, the type carries certain p-adic discrete valuations 

Vi ■ Qpix)* —> 1, and semigroup homomorphisms, 

w,: Zp[x] \ {0} ^ Z>o, P{x) ^ ord^,_,(i?,_i(P)), 

where RQ{P){y) G Fp[y] is by definition the reduction modulo p of P{y) /p^'^^\ 
For instance, the integers wi(/), . . . are the above mentioned exponents 

ni,..., Ur+l. 

To avoid confusion, in case of working simultaneously with different types, we 
add a superscript with the type to every component or datum: A*, e\, etc. 

Definition 2.1. A type t of order r is said to he complete if ujr+i{f) = 1- 

Theorem 2.2 f (HN[ Sect. 3]). Lett be a type of orderr such thatur+i := ujr+i{f) > 
0. Then, t singles out a monic p-adic factor ft{x) G Zp[a;] of f{x), uniquely deter- 
mined by 

nt := deg/t = TO^+in^+i, ujr+iift) = "-r+i- 
If t is complete then ft (x) is irreducible in Zp [x] , so that t actually singles out 
a prime ideal p dividing p TLk ■ The ramification index and residual degree of p are 
given by: 

e(p/p) = ei • • • e^, /(p/p) = /o/i • • • /r- 

Therefore, the branches whose type t is complete are end branches that bear 
no further ramification. Whenever the algorithm constructs a type, it is sent to 
a list stack, if it is non-complete, or to a list completetypes if it is complete. 
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The algorithm runs till the list stack is empty, or equivalently, till all types are 
complete. 

Along the analysis of the branch labelled by t, the algorithm computes the 
positive integer 

indt(/) :=/o---/.-iind(7V-(/)), 

where iiid{N,^ ( f)) counts the number of points with integral coordinates that lie 
below or on the polygon, strictly above the horizontal line that passes through the 
last point of N^{f) and strictly beyond the vertical axis. 

Definition 2.3. For any natural number r > 1, we define 

ind,(/) i^dt(/), 
tet,_i(/) 

where tr-i (/) is the set of all types of order r — 1 that are computed by the algorithm. 

Theorem 2.4 f [HNi Sect. 4]). There is some r > such that all types oftr{f) are 
complete. In this case, 

(2.1) ind(/) = indi(/) + ...+ind,(/). 

In particular, by a simple accumulation of all indt(/) to a variable totalindex 
the algorithm computes ind(/) as a by-product. 

Finally, the algorithm incorpores a crucial optimization. A refinement process 
looks for an optimal choice of the representative (f>r+i{x) of the type t, for each of 
the ramified branches of order r + 1 that sprout from the branch represented by t. 
This process requires the introduction of a variable H — _ff * called the cutting slope, 
attached to each type. Initially one takes H = 0, and, if (pr+iix) is not optimal 
for certain subbranch, then for the type t' that will be responsible for the analysis 
of the subbranch, iJ* is given certain positive value. This refinement causes a 
slight modification in the way the higher order indices indt(/) are accumulated to 
totalindex 

The main loop analyzes a type from stack and determines its ramification. 

Main loop of Montes' algorithm. At the input of a type t of order r — 1, for 
which LOrif) > 0, with representative (f>rix) and cutting slope H: 

1. Compute the partial polygon, N^{f), that gathers all sides of slope less than 
—H of the Newton polygon of 7--th order Nr{f), and accumulate to totalindex 
the value 

indf (/) /o • • ■ fr-i (^ind(7Vf (/)) ~ ^Hiil - 1)) , 
where £ is the length of N^{f). 

2. FOR every side of N^{f), let A < -H be its slope, and do 

3. Compute the residual polynomial of r-th order, R\{f){y) G ^r[y], and fac- 
torize this polynomial in [y] . 

4. FOR every irreducible factor ip{y) do 

5. Compute a representative (pr+iix) of the type t' :=(... , 4>r{x); X, ip{y))- 

6. If t' is complete, then add t' to completetypes, and continue to the next 
factor of Rx{f)iy). 
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7. If dcgif,' = 1 and A G Z (the type must be refined), set <~ <j)r+i{x), 
H*' ^ A, add t to stack and continue to the next factor of R\{f){y). 

8. (Build a higher order type) Set ^ 0, add t' to stack and continue to 
the next factor of Rx{f){y). 

After each iteration of the main loop, the set T = stackU completetypes, of all 
types produced by the algorithm determines a factorization of f{x) in and 
the corresponding partition of the set V of prime ideals of K lying above p: 

m = n ^ = II ^t, 

ter teT 

where Vt denotes the set of prime ideals of V such that the corresponding p-adic 
irreducible factor of f{x) divides ft{x) in Zp[x]. Clearly, 

nt := deg/t = ^ e(p/p)/(p/p), n = ^ nt- 
peVt ter 

3. Quotients of 0-adic developments 

Definition 3.1. Let 4>{x) G Z,[x] be a monic polynomial, and let 

}{x) = ao{x) + ai{x)4>{x) H h as{x)4>{xY , 

with ai{x) e Z,[x], degai{x) < deg(l){x), be the (j)-adic development of f{x). 

The quotients attached to this 0-development are, by definition, the different 
quotients qi{x), . . . , qs{x) that are obtained along the computation of the coefficients 
of the development: 

f{x) = (j){x)qi{x) + ao{x), 
qi{x) = (f){x)q2{x) + ai{x), 

qs^i{x) = (t>{x)qs{x) + as^i{x), 
qs{x) = (l){x) ■ + as{x) ^ as{x). 

Equivalently, qi{x) is the quotient of the division of f{x) by (p{xy; we denote by 
ri(x) the residue of this division. Thus, for all 1 < i < s we have, 

(3.1) f{x) ^ ri{x) + qi{x)(l){x)\ ri{x) ^ ao{x) + ai{x)4>{x)^ t-ai_i(x)0(x)'"^ 

The aim of this section is to compute, for certain quotients q{x) attached to 
certain 0-adic developments of /(x), the highest power of p such that q{6)/p"^ 
is integral. To this end we shall need three lemmas. The first one is an easy 
generalization of |HN[ Prop. 2. 7], in a spirit similar to that of [HNl Lem.4.21]. 

Lemma 3.2. Let t = {(f)i{x); ■ ■ ■ , 4)r{x)', Xnipriu)) be a type of order r . Let 

J = {(ii, ...,>)£ I < ji < e,/„ 1 < ^ < r}, 

and for any j G J denote <i?(x)-' — (/)i(x)^^ • • • (jfrixy^' . Let P{x) G be a polyno- 
mial of degree less than m^+i, and consider its ^-multiadic development: 



P{x) = degflj < mi, Vj G J. 

Then, Vr+i{P) = mmiej{vr+i{aj{x)^{xy)} . 
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Proof. Since Vr+i is a valuation, it is sufficient to check that Vr+i{aj{x)^(xy) > 
Vr+i{P), ioT all j G J. Let us prove this inequality by induction on r > 1. For 
r = 1 this is proven in [HNi Prop. 2. 7]. Suppose that r > 2 and the lemma is true 
in order 1, . . . , r — 1. Consider the (/)r-adic development of P{x): 

0<i<e,./. 

By (HNl Prop.2.7], 

(3.2) Vr+l{P) < Vr+l{Pjix)(f)r{xy) = CrVriPj) + jVr+l{(j>r), 

for all < j < Crfr- Consider now the set 

J° = {{Ji, ■ ■ • e N''"i I < ,n < e,/„ 1 < z < r}, 

and for any < j < Crfr consider the embedding 

J° J, j° = (ji, . . . , ^ (j°, j) := (ji, . . . , jr-l, j)- 

The (^1, . . . , 9ir-i)-multiadic development of Pj{x) is: 

and by induction hypothesis: 

Vr{Pj) <Vr{aQO^j){x)(f>i{xy^ • • • (/(^-l (x)^""' ) 

= t;r(a(jOj)) +ilWr(0l) H 1- jr-lWr(0r-l), 

for aU e jO. By [HN, Prop.2.7,(l)] and [HNl Lem.2.2,(2)] we have 

erWr(a(jO,i)) = ^^r+l(a(jOj)); erVr{(j)i) = Vr+l{4'i), VI < 1 < T, 

so that the last inequality, together with (|3.2|) gives the desired inequality: 

Wr+l(-P) < «r+l(a(jO,j)) + jlWr+l(</'l) H h 1 Wr+1 ((/"r- 1 ) + jVr+l{4'r), 

for aU j" e J" and all < j < e^/^. □ 

We shall also use the following elementary combinatorial result. 

Lemma 3.3. Let Mi, . . . , Mr be positive integers. Then, the set 

{ji + kMi + ■■■+ >MiM2 • • • Mr-i I e Z n [0, Af,), 1 < i < r}, 

coincides with the set Z n [ 0, M1M2 ■ ■ ■ Mr). □ 

For any pair of indices, s, t, such that s < t < r, let us denote 

As^t ■= Gs+i • • • er-i{esfs) ■ ■ ■ (et_i/t_i). 

With this notation, the closed formulas for Vr{<t>t) and Ws(0s) from |HN| Prop. 2. 15] 
can be written as: 

«r(0t) = Y!k=x Ak,thk, t <r, 

(3-3) Wr(0r) = Yl^l ^k,rhk, 

esAs,tVs{(f>s) = Z]fc=i Ak,thk, s <t <r. 
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Lemma 3.4. Suppose that s < r, and let {js, ■ ■ ■ , jr-i) G N'' ^ be such that < 
jt < etft, for all s < t < r. Then, 

s<t<r 

Proof. The inequality is equivalent to 

X! jt{esfs) ■ ■ ■ (et_i/t„i) < (Cs/s) • • • (er-lfr-l), 
s<t<r 

and this is an immediate consequence of the previous lemma. □ 

Corollary 3.5. Let r > 2, and let t = - ,(t>r-i{x); Xr-i,^r-i{y)), be a 

type of order r — 1, with representative (j)r{x). Let (ji, . . . , jV-i) G be such 

that < jt < Ctft, for all 1 < t < r. Then, 

jlVr{(f>l) H \- jr~lVr{(l>r-l) < I'rl^r)- 

Proof. By p.3p . we want to show that 

t r-l 
l<t<r fe=l fe=l 

For each value of 1 < fc < r, the terms involving hk in both sides satisfy this 
inequality by the previous lemma: 

itAk^hk < Ak^rhk. 

k<t<r 

□ 

Proposition 3.6. Let r > 1 and let t — (0i(a;);-- - ,(f)r~i{x)\\r~i,'ipr-i{y)), be 
a type of order r — 1, with representative (f>r{x), such that Ur :— i^r{f) > 0. Let 
1 < i < Ur be an integral abscissa, and let Yi ^ Q be the ordinate of the point of 
(f) with abscissa i. Let qi{x) be the i-th quotient attached to the (pr- development 
of f{x). Then, for every prime ideal p of K lying above p we have: 

e(p/p) ~ eoei • • -e^-i ' 
where cq = \ by convention, and fp is the discrete valuation of K determined by p. 

In particular, if Hi :— \{Yi — ivriiir)) I {eoCi ■ ■ ■ er-i)\, then qi{6)/p^' belongs to 
Zk. 

Proof. Let p € V he a, prime ideal of K lying above p, and let Fp{x) S Zp[x] be 
the monic irreducible factor of f{x) corresponding to p. Let us choose a root 6"^ 
of Fp {x) ; consider the embedding Lp-. K ^ Qp determined by sending to 6"' . If 
we denote, Wp := e{p/p)~^Vp, then we have: Wp(a) ~ v{Lp{a)), for all a G K. In 
particular, for any polynomial P{x) G Z[x], 

(3.5) Wp{P{e))=v{P{9n). 

Let Ao be the slope of the side S of N^{f ), whose projection to the horizontal 
axis contains the abscissa i. If i is the abscissa of a vertex of N~{f), then we let 
S be one of the two adjacent sides. The prime ideal p satisfies one and only one of 
the following conditions: 

(1) peVf, fort' = - , A, V(2/)). 
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(2) p e Vf, for t' = {<j)i{x); ■ ■ ■ ,(t>s{x); \,^p{y)), for some 1 < s < r, with 

fp{y) i^s{y) if A = As. 

(3) p\MO) inZK. 

We shall prove the inequality (j3.4p by an independent argument in each case. 
In the two first cases, since Fp{x) is a factor of ft'{x), the element 6"' will be a 
root of ft'{x) too. This allows us to apply the results of [HNl Sect. 3] to the pair 
ft'{x), 6"'. We shall denote throughout the proof: e = eoCi • ■ • e^^i and vt = , for 
all 1 < i < r. 

Case 1: peVf, for t' = (0i(x);--- , 0r(x); A, ^(y)) 

By p.ip . the Newton polygon N~{f ) splits essentially in two parts: N~{ri) and 




This is not always true because, depending on the values of Vr{ai~i(t>^~^) and 
Vriaicf)"^), the two parts of the side S in the polygons N~{ri) and N^{qi(f)\.) might 
change. The next figure shows different possibilities for these parts of S. However, 
the line Lx^ of slope Aq that first touches both polygons from below is still the line 
determined by S. 
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If |A| > |Ao|, we apply [HN| Prop. 3. 5, (3)] directly to the polynomial qi{x)(j)r{xy 
to get: 



Y^ + i\X\ 



If |A| < |Ao|, we apply I HN l Prop. 3. 5, (3)] to the polynomial ri{x) and we get the 
same result, thanks to (|3.ip : 



Y, + i\X 



Yi+i\Xo\ 




\M > |Ao| 



By the theorem of the polygon [HNl Thm.3.1], and (^3)1 : 

w^iMd)) = viMO^)) - + |A|)/e, 

and we get the desired inequality: 



Y, + . Vr{(f>r) + |A| _Y,~ iVr{(f>r) 



Case 2: p e Pf, for t' = ((/)i(x); • • ■ , (psix); X,ip{y)), for some 1 < s < r, with 

V'(y)^V's(y) if A = As 

Recall that ri{x) = J2o<j<i^j(^)'^rixy ■ Our aim is to show that Wp{aj{0y is 
sufficiently large, for all < j < i. More precisely, we want to show the following 
two inequalities: 



(3.6) 



Wp{(t>r{0y < 



Vr{(j}r) 



(3.7) ^^[a,{e))>^-^-(^-^-w,{ct>r{ey], VO<j<z. 



Let us show first that these two inequalities imply the desired inequality (j3.4l 
In fact, from ()3.6|) and ()3.7|) we get: 



(3.8) yJ^[a,{ey>'^-{^-J)(^-^^w,{<t>r{e))], yo <j<i, 
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and we can argue as follows: 



> min {wp{aj{e)(l3ri0y)} - iwp(<?!)r(6')) 

0<j<i 



> 



. ^Vr{(f>r)_ ^ Yjg - iVr{4'r) ^ Yj - iVr{4>r) 



Let us compute now the exact value of Wp{(j)t{9)) — v((/)t(6'f )), for all s < t < r. 
For t = s, the theorem of the polygon [HN[ Thm.3.1] shows that: 



(3.9) 



Vs{^s) + \X\ 



eoci • • ■ Cs-i 

For s < t < r, v{(f>t{0^)) can be computed by applying [HN[ Prop. 3. 5, (3)] to the 
polynomial <j)t{x). We know that Ns{(j)t) is one-sided with slope A^; thus, there 
are two different situations according to |A| > |As|, or |A| < |As|, reflected in the 
following figures: 



Ns{ct>t) 



Vs(4>t) 




Ns(<l>t) 



Vs{4>t) 




If A ^ As, the A-component of Ns{(t>t) reduces to a point, and the residual 
polynomial of s-th order R\{(f>t){y) is a constant [HNl Defs.1.5,2.21]. If A = As, 
then R\{(j)t){y) is a power of ?/'s(y); since ^ i^s{y), we see that R\[cf)t){y) is not 
divisible by 'ip{y) either. Thus, [HN, Prop. 3. 5, (5)] shows that v{<j)t{0^)) is always 
equal to H / {e^ei ■ ■ ■ es_i), where H is the ordinate of the point of intersection of 
the vertical axis with the line of slope A that first touches Ns{4>t) from below. By 
taking a glance at the figures above we see that: 

(3.10) wp{(l)t{0)) = v{(l>t{9^)) ^ , s<t<r, 

eoei • ■ • es^i 

where ^ := min{|A|, |As|}. Since A^s(0t) is one-sided with last point {mt/ms, Vs{4>t)) 
[HNl Lem.2. 17,(1)], we have 

v,{^t) = «s(C/'"0 = (mt/ms>s(0s) = (es/s) • ■ • (et_i/t_i)«s(0s). 
so that (j3.10p can be rewritten as 

(3.11) v{UO')) ^ {esfs) ■ ■ ■ (et-i/t-i) ""'^^'^ ^ ^ , s < t < r. 

eoei ■ • ■ es_i 
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Since |A| > /i, we deduce from dHJ]), (j3lT|) and ^3^ : 

min{i,r-- 1} 

(3.12) ev{(t)t{6''')) > esAs,t{vs{4's) + fJ-) ^ Vr{(f)t) - ^ Ak,thk + SsAs^tfJ', 

k=s 

for all s < t < r. 

Let Jg :— {{js, ■ ■ ■ ,jr-i) & I < jt < etft, s < t < r}, and consider the 

following multiadic development of aj{x): 

%(^) = X! ki^)'l^'^i^y' ■ ■ ■ (l)r-i{xy'"\ degfcj < nis. 
Fix a multiindex j G Js such that 

^v{bj{9^)) + jMM0')) + ■ ■ • + >-i«(0.-i(0P)). 

By [HNl Prop.2.9], w(6j(6'P)) = 'yr(6j)/e. Therefore, by KWi we get 
ei«p(a,((?)) = e«(a,(0P)) 

s<i<r V A;=s+1 / 

> Vr{aj) + ^ jf L4s,f(es/^-/is)- ^ ^fc,t/ifc I • 

s<i<r \ fc=s+l / 

the last inequality by Lemma On the other hand, for t = r we have by (|3.12p : 

r-l t 
Vr{(l)r) - eWp{(l)r{9)) = -esAs,rfl + ^^Ak,rhk = As,t{hs - esU) + ^ Ak,thk>0, 

k—s k—s+1 

because hg > Cg/i. This proves (|3.6p . Also, (|3.7[) is reduced now to: 



X! Jt Ms,t(/is-es^)+ ^ ^fe,t/ife I <^Ak,rhk - e^As^rl^- 

s<t<r \ k=s+l / k=s 

If we consider only the summands involving hs and ^, we get the desired inequality 
by Lemma [33] 

{hs - Cs^^) ^ jtAs,t < {hs - es^J,)As,r■ 

s<t<r 

Finally, if we consider the summands involving hk, for any fixed s < k < r, we get 
again the desired inequality by Lemma 13.41 

hk ^ jtAk,t < hkAk^r- 

k<t<r 

Case 3: p \ <l)i{0) in Zk 

We have now Wp{(t>t{9)) — 0, for all 1 < t < r. Thus, (|3.6p is obvious, and 
arguing as in Case 2, it is sufficient to check the inequality (|3.7p . Let us consider 
the multiadic development of aj{x): 

^ji^) = '^hi^)^ii^y^ ■ ■ ■ 4'T-i{xy''-^ , degfoj < mi. 
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We fix j e J such that v(aj{9P)) > v{bi{9^)). By jHH Prop.2.7,(l)], Vr{bj) = 
evi{bj), so that 

ewp{aj{9)) ^ev{aj{0^)) > et;(6j(6'P)) > evi{bj) = Vrib^) > (by LcmmaE^]) 

>Vr{aj) - jlVr{(f>l) jr-^lVr{(f>r~l) > Vr{aj) -Vr{(j>r). 

the last inequality by Corollary [331 D 



4. Computation of a p-integral basis 

4.1. Computation of n integral elements. Let t = {(j)i{x); ■ ■ ■ ,4>r{x);Xr, i^r{y)) 
be a type of order r. We recall that the reduced type determined by t is by definition 

t° {(f>i{x);- ■ ■ ,(j)rix);Xr). 

Note that t° is not actually a type. 

Definition 4.1. Let t° — (0i(a;); • • • , (j)r{x); Xr) be a reduced type of order r. We 
say that is regular (or more properly: /-regular^, if the r-th order residual 
polynomial R\^{f){y) G is separable. 

Imagine that in one iteration of the main loop of Montes' algorithm, we analyze 
a type t of order r — 1, with representative 4>r{x). We consider one of the sides S 
of N~ (f), with slope (say) A,-, and we factorize the r-th order residual polynomial 
in ¥r[y]: 

(4.1) RxAf)iy)^^rAyr---^rAyr'. 

Then, we are about to ramify t into g types of order r 

U ■■= {(f>l{x); ■■■ , (j)r{x);Xr,^r,i{y))^ ^ <i< 9, 

except for the branches that might be refined (if G dcgipr,i — 1 and t^ is not 
complete) . 

Remark 4.2. If the reduced type (</>i(a;); • • • ,(j>r{x);Xr) is regular, then ti, . . . ,tg 
are all complete. Moreover, in this case 

nti H h «tj, = "ii(ei/i) ■ • • {er-ifr-i)iiS), 

where £{S) is the length of the projection of S to the horizontal axis. 

Suppose that CT completetypes is the output list of complete types provided 
by Montes' algorithm, and let CT''''^ be the set of all regular reduced types t°, for 
t S CT. Let us show how these two sets can be used to construct a p-integral basis. 

We introduce some notation attached to any t = {(j>i{x); • • • , (jfrix); Xr,ipriy)) S 
CT of order (say) r. For all 1 < t < r, let [at, bt] be the projection to the horizontal 
axis of the side of slope At of Nf{f). Now, for all integer abscissas, at < j < bt, 
we denote by Yt j the ordinate of the point of Nf{f) of abscissa j, and 
(4.2) 

1 , if t < r and j = bt, 

j-th quotient of the 0t- development of f{x), otherwise. 



Qtjix) 



3, if < < r and j ^ bt, 

(4.3) Ht^.ix) = { Yt,j - ]vt{4>t) 



otherwise. 



eoei • • • et-i 
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By Proposition [SH qt,j{9)/p^"'-' ^ is integral, for all t,j. We shall construct a set B 
of n elements of Zk ; the basic idea is that each t G CT contributes with nt elements 
to the basis B. However, when for some t e CT, the reduced type t'^ is regular, 
then we compute the contribution of the complete types ti, . . . ,tg of Remark 14.21 
altogether. 

Case t° regular. We take all complete types ti, . . . , € CT, that have a common 
reduced type, t'^ = ■ • ■ , (firix); A, ), which is regular. Let J be the set: 

{(jo, ■•■,>)€ N'+i I < jo < mi; h ~ ctft < jt <bt, I <t < r; at < jr < K}. 
In this case, we add to B the set: 



plHr,j,.^ ^-f^lJl J 



(jo, . . . , jr) e J 



Note that | Qt" I — "T-i(ei/i) • • • (er-i/r-i)(fcr — Or) = "ti + • • ■+nt^ , the last equality 
by Remark 1321 

Case t*^ not regular. We take a complete type t € CT, such that t° is not 
regular, and we compute a representative ^^+1(2;) of t. The polygon N~^j^{f) has 
length one, with integer slope A^+i, and its last point is (1, ■i;r+i(/))- We compute 
the first quotient of the 0r+i-adic development of f{x): 

f{x) = (l)r+i{x)q{x) + a{x), dega(x) < m^+i. 

By Proposition 13. 6[ if H :— (wr+i(/) — Vr+i{(t)r+i) / {eoCi ■ ■ ■ Cr) , then q{6)/p^^^ 
belongs to 'Zk- We add then to B the set: 

'q{9)qr,jM---qi.ni0)O'" 



Qt ■■= 



jt e {bt - etft,bt], I <t<r; Jq e [0,mi) 



Caution: qt,j{x) are now considered as quotients of ^(-developments with respect 
to a type of order r + 1; thus, following the criterion of (|4.2p and l|4.3p . we take 

qr,br.{x) = 1, i?r,&, = 0. 

Note that |Qt| = "ii(ei/i) • • • (e^/r) = nt. 
Conjecture 4.3. The following family is a p-integral basis of K: 

B={ U 2to)u( U ^ 

4.2. Montes' algorithm suitably modified to compute a p-integral basis. 

Montes' algorithm computes already the quotients of the different (/)t-adic develop- 
ments of f{x) along the process of determining the coefficients of the development. 
Thus, the algorithm can be easily adapted to compute a p-integral basis, just by 
computing the exponents Ht,j of the denominators and the adequate products of 
quotients to build the families Qt" and/or Qf 

A part of these products is shared by the subbranches of a type, when it ramifies; 
therefore it is more efficient to compute the partial products before a type ramifies, 
so that all subbranches can share these partial products without repeating their 
computation. To this end, we attach to each type t two lists ProductOf Quotients 
and ProductOf QuotientsValuations, denoted PQ = PQ*, PQVals = PQVals* in 
the sequel for simplicity. These lists are initialyzed by 

PQ* = {l,a;,...,a;'^'=s0i-i|^ PQVals* = {0, • ■ • ,0}, 
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and they are enlarged in the foUowing way. Suppose that we analyze a type of 
order r — 1 from the list stack: t = (0i(a;);--- , (j)r-i{x); Xr-i,tpr-iiy)), with 
representative (t>rix)', we compute the Newton polygon N~{f), we consider one of 
its sides S with slope (say) = —hr/cr, and we factorize the residual polynomial 
as in (|4.ip . If the reduced type (01 (x); • • ■ , (j)r-i{x); Xr-i, 4>r{x); \r) is not regular, 
we compute fmax :— maxi<i<g{deg V'r.i} and we construct the following two lists: 

(4.4) EnlargedPQ ^ PQ* U qr,b,-i ■ PQ* (mod f{x)) U • • • 

U <?r,6,-e,fmax • PQ* (mod/(x)), 

(4.5) EnlargedPQVals ^ PQVals* U + PQVals* U ••• 

J. fmax + PQVals'^ . 

For eac branch t' — {{(f>i{x); • ■ ■ ,(l)r-i(x)]Xr-i,4'r(x)-,Xr,'4>r,i{]j)), that sprouts 

t' t' t 

from t, the lists PQ , PQVals are the sublists of the first e^/rjiPQ elements, 

respectively of EnlargedPQ, EnlargedPQVals. 

Finally, we introduce two global lists numerators, denominators. When the 

algorithms ends, the p-integral basis is made up by all elements qi{0) j^f ^ , 1 < i < f^, 

with qi{x) G numerators and vi G denominators. The main loop of Montes' 

algorithm is then modified in the following way: 

Main loop of Montes' algorithm, modified to compute a p-integral basis. 

At the input of a type t of order r — 1, for which uJr{f) > 0, with representative 
4>r{x) and cutting slope H: 

1. Compute the partial polygon, N^{f), that gathers all sides of slope less than 
—H of the Newton polygon of r-th order Nr{f), and accumulate to totalindex 
the value 

indf (/) := /o • • • fr-i (^ind(7Vf (/)) - ^Hi{e - 1)^ , 
where i is the length of N^{f). 

2. FOR every side of N^{f), let Xr < ~H be its slope, and do 

3. Compute the residual polynomial of r-th order, R\^{f){y) G IP'rly], and 
factorize this polynomial in Fr[2/]. 

3.1 If Rxr.{f){y) is separable, then 

3.1.1 Add <j<b 1r,jPQ (mod /(x)) to numerators. 
Add <j<b l-^r.j + PQValsJ to denominators. 

3.1.2 Add the types t' = {(j)i{x);--- ^ (t)r{x); Xrjipriv)) to completetypes, 
where ijjriy) runs on all irreducible factors of R\^{f){y). Continue to the next side 
of7V/^(/). 

3.2 Build the lists EnlargedPQ, EnlargedPQVals, as indicated in ((44)) and (|43)) . 

4. FOR every irreducible factor ?Ar(y) do 

5. Compute a representative (j)r+i{x) of the type t' :=(... , (j)r{x)] Xrjijjriy))- 

6. If t' is complete, then add t' to completetypes, and 

6.1 PQ ^ {EnlargedPQ(i) | 1 < i < e^/J, 

PQVals ^ {EnlargedPQVals(i) | 1 < i < Crfr}- 
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6.2 Compute the first quotient q{x) of the (/)r+i-adic development of 
f{x), and the exponent H given by Proposition 13.61 Add q{x)PQ (mod f{x)) to 
numerators and + PQValsJ to denominators, and continue to the next factor 

ofi?A(/)(y). 

7. If degip = 1 and A e Z (the type must be refined), set (f>l{x) ^ (j)r+i{x), 
^ A, add t to stack and continue to the next factor of Rx{f){y). 

8. (Build a higher order type) Set i?*' ^ 0. 

8.1 PQ ^ {EnlargedPQ(i) \l<i< e^/^}, 
PQVals ^ {EnlargedPQVals(i) | 1 < i < e^/^}. 

8.2 Add t' to stack, and continue to the next factor of R\{f)(y). 



When the algorithm ends, the list numerators contains n polynomials qi{x) with 
integer coefficients, and the list denominators contains n non-negative integers Vi. 
We explain now how can be checked that the elements qi{0)/p'^^ form a p-integral 
basis indeed. 

Consider the Z-modules 

(g^(6l) I q,{x) e numerators)^ C Z[9], 

O := ^^-77-^ Qiix) G numerators, iyi{x) € denominators^ C Z^, 

We compute ind""™ v {Z[d]: 0""™) as the p-adic value of the determinant 
of the matrix made up with the coefficients of all qi{x). Also, we can compute 

v{0: O"™) =Er=i'^- 

On the other hand, since ind(/) — v{'Lk'- ^[6*]) is computed by the algorithm, 
we know the value of v {Zk : C"""") ind"""' + ind(/). Thus, B is a p-integral basis 
if and only if 

n 

(4.6) ^i^, =ind"™ + ind(/). 

i=l 



4.3. Triangular basis. We say that a p-integral basis is triangular if it consists 
of n integral elements of the type qi{9)/p'^% with qi{x) € Z[x] monic of degree i. 

It is well-known that the Chinese remainder theorem can be used to construct an 
integral basis of Zk from triangular p-integral bases for the different prime divisors 
p of disc(/). 

Unfortunately, the p-integral basis B = {qi{0)/p'^' | 1 < i < n}, produced by the 
above modification of Montes' algorithm, is not triangular. On the other hand, for 
n large the computation of the determinant of the matrix of the coefficients of the 
numerators qi{x) is inefficient. We can save these two handicaps by computing a 
p-stem of a p-integral basis from the output list B. We recall that a p-stem is a 
family of generators of Zk ®z Z(p) as a Z(p) [6']-module, of the form qi{9)/p^\ with 
monic polynomials qi{x) G Z[x\ having minimal degrees. 

The p-stem is obtained by a single call to the MAGMA routine GroebnerBasis 
in order to compute Groebner generators of the ideal of Z[x] generated by {p^} U 
Sjf-'^iq^i^x) I 1 < i < n}, where v — maxiji/i}. The output of this routine is a list of 
polynomials p"^, p'^^'^'^ gi{x), . . . ,p'^^^^'' gk{x), with gi{x) £ Z[x] monic polynomials 
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having strictly increasing degrees; the p-stcm of I^k is given by 

, 51W ... 

and this is the final output of the algorithm. If di := deggi. for 1 < i < fc, the 
triangular p-'mtegral basis determined by the p-stem is: 



pf^l p^l pf^k pl^k 

Once the p-stem is computed, in order to show that we obtained a p-integral 
basis it is sufficient to check that 

(4.7) {d2 - di)iii + • • • + (n - dk)nk = ind(/). 

If this equahty holds, the algorithm prints a message of maximal order, and other- 
wise a message of non-maximal order. 

This results into a more efficient computation, but the routine GroebnerBasis 
uses much more memory. 



4.4. An example. Let us show how the algorithm works with an example. Take 

p = 2 and f{x) = + 4,7.6 ^ iq^3 

Since f{x) = x^"^ (mod 2), the set to(/) contains only one type to of order 
zero, represented by the irreducible polynomial ipo{y) — y; we take 4>i{x) = a; as 
a representative of this type. The Newton polygon of first order of f{x) has two 
sides, with slopes —2/3 and —1/3, and indi(/) = ind(iVi(/)) = 23. 




The residual polynomials are: 

R-2/3{.f)iy) ^y^ + y + l, R-i/sifM = iy + i)'- 

Thus, the type to ramifies into two types of order one: 

t = {x; -2/3, y^ + y + 1), t' = (x; -1/3, y + 1), 

and ti(/) = {t.t'}. The type t is complete, and it singles out a prime ideal pi 
with e(pi/2) = 3, /(pi/2) = 2. Since its reduced type t° is regular, the six first 
quotients of the (?ii-adic development of f{x) determine six elements from the final 
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list B: 



nix) 


= +4a;'5 + 16a;2, 


Hi 


= 16/3, 




gi(0)/25 


Q2{x) 


= + 4a;'* + 16a;, 


H2 


= 14/3, 




Qsix) 


= x^ + 4x3 + 16, 


H3 


= 4, 


•s^ 


q4{x) 


= X* + 4x2, 


Hi 


= 10/3, 


-w 


qd0)/2^ 
qd0)/2^ 




= x'' + 4x, 


H5 


= 8/3, 


-w 




= x6 + 4, 


He 


= 2 


-w 



The type t' is not complete, and it requires to do some more work in order two. 
Before analyzing the future ramification of t' we enlarge the list PQ* with the 
quotients qwix) — x^ and qii{x) = x of the 0i-adic development of /(x): 

PQ*' ^ {l,x,x2}, PQVals*' ^ {0,1/3,2/3}. 

Let us choose (j)2{x) = x'^ + 6 as a representative of t'. The (/)2-adic development of 
fix) is: 

fix) = cl)2ix)'^ ~ 2402(a;)3 + 220(/)2(x)2 - 89602(x) + 1408. 

The principal part of the Newton polygon of second order has length 2 = cj* (/); 
hence we have to draw only the first three points to determine N^if)- Since 
1^2(02) = W2(2) = 3, we have 1)2(1408) = 21, W2(89602) = 24, ?;2(22O0i) = 12, so 
that N2 if) has only one side of slope —9/2, and ind2(/) = ind(A'^2 (/)) = 4: 



21 1 












12 ■ 
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The residual polynomial of second order is already separable: i?_9/2(/)(j/) = 
y + 1. Thus, t' is extended to a unique type of order two: 

t" = (x; -1/3, x3 + 6; -9/2, y + 1), 

which is already complete. It singles out a prime ideal p2 with e(p2/2) = 6, 
/(P2/2) = 1. Since the reduced type (t")° is regular, the two first quotients of 
the (/)2-adic development of /(x) determine six elements of the final list B. The two 
quotients and the exponents of their denominators, given by Proposition [321 are: 

Q^{x) = ^2(a;)3 - 24^2(2;)^ + 220^2(2;) - 896 = x^ - Qx^ + 40x3 - 224, 

Hi = (16.5 -3)/3 = 9/2, 

Q2(x) = <f)2ix)'^ - 2402(a;) + 220 = x^ - 12x3 + 112, H2 = (12 - 6)/3 = 2. 



HIGHER NEWTON POLYGONS AND INTEGRAL BASES 



19 



By multiplying them by PQ* , each one determines three elements of B: 
Qi{e) Qi{0)O Qme^ Q2{9) Q2{0)9 Q, 



24 ' 24 ' 25 ' 22 ' 22 ' 22 ' 
and our list B is complete, with 12 integral elements. 

The output list of complete types is CT = {t,t"}, parameterizing the prime 
ideals pi, p2- The index ind(/) = v^Ek ■ '^[0]) is equal to: ind(/) = indi(/) + 
ind2(/) = 23 + 4 = 27. Since t;(disc(/)) = 69, we get ?;(disc(i4:)) = 15. 

The reader can check that the sum of the exponents of the denominators of 
the elements in B is 39, and ind""™ = 12; thus, 6 is a p-integral basis indeed, by 
criterion (|4.6p . 

However, the algorithm does not check this. It outputs directly the p-stem: 
6*3 6i5 + 202 q6 ^ 26»3 + 26*5 + 86*2 6*" + 46^ + 169^ 
Y' 22 ' 23 ' 24 ' 25 ' 

by a call to the GroebnerBasis routine, as explained above. Then, it checks that 
criterion (14.71) is satisfied: 



2-l + l-2 + 2-3 + 3-4+l-5 = 27; 
thus, it outputs a message of maximal order. 

4.5. Proof of the conjecture when f{x) is irreducible over Zp[a;]. In [EN09] 
it is proved that S is a p-integral basis if all complete types of the output list CT 
have order less than or equal to one. We deal now with the case f{x) irreducible 
over Zp[x], or equivalently, the output list CT consists of only one type. 

Theorem 4.4. Supose that f(x) is irreducible over Xplx]. Let 

t = {(j)i{x); ■ ■ ■ ,(l)r{x); A,.,-0r(y)), 

be a complete type of order r, and denote Ui := LOi{f), for 1 < i < r+ 1. Let Cq = 1, 
fa = mi and J := {(jo, .-.,>) G N'^^^ | < < e^/^, for all < i < r}. Then, 
with the notation of ^-2^ and ^-3^, the family of all 



(A fi] Qr,nr~jA9) ■ ■ ■ qi..m-ni9)9^° , . J 

is a p-integral basis of K . 

Proof. By the Theorem of the polygon and the Theorem of the residual polynomial 
[HNl Thms.3.1,3.7], for each l<i<r: 

(1) Ni{f) is one-sided, with slope and length = {cifi) ■ ■ ■ (crfr), 

(2) ind,(/) =/o---/.-iind(7V,(/)), 

(3) R,{f){y) ^ Uvr-^' inF,[y]. 
Since t is complete, we have n^+i = 1, and 

n = miUi = nt = /o(ei/i) • • • (e^/r), 

for all 1 < i < r + 1. Thus, the family of (|4.8p consists of n algebraic numbers, all 
of them integral by Proposition 13.61 Let O C Zk be the sub-Z- module generated 
by these elements. 

By (|3.ip . the degree of each qi^m-jiix) is equal to n — {ui — ji)mi = jirrii. By 
Lemma 13. 3| the degrees of the numerators of the n integral elements of (|4.8p take 
the values 0, 1, . . . , n — 1; therefore, the elements qr,n^-j^ {9) ■ ■ ■ 9i,ni-ji {9)9^° , for 
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(jo, ■ ■ ■ ,jr) e J, are a Z-basis of the order Z[9]. In particular, O is a free Z-module 
of rank n that contains and 

V {O: lie]) = lHr,n^-Jr +■■■+ Hl,rn-n\ 

(jo,---,ir)eJ 

= fo X! [Hr^n,^-],^ ^ |-^l,ni-JiJ- 

{o,ju...,jr)eJ 

Therefore, in order to show that v{Zk : O) = it is sufficient to check that 
v{0: 1[9]) = v{Zk: ^[O]) = ind(/). By this is equivalent to 



(4.9) /o VHr,n.-i,. 
{o,3i,...,jr)eJ 



■ ■ + i/i,„i-jj = indi(/) + • • • + indr(/). 
By a simple glance at the Newton polygon Ni{f): 




we see that Vi{f) = (</»"') and = v^{f) + {ui - j)\Xi\; hence, 

Yi.7H-j, - {"n-i - ji)vii4>i) _ Viif) - {rii - ji)vi{(f>i) 



eo • • ■ ei_i 



eo ■ ■ • ei-i 



eo ■ • ■ ei-i 



the last equality by the Theorem of the polygon [HN[ Thm.3.1]. 
Therefore, the equality (|4.9p can be rewritten as 



(Oji,...,>)e./ i=i 



indr(/), 



and this equality is proved in |HN[ p. 55], as one of the final steps of the proof of 
the Theorem of the index. □ 



5. Some numerical examples 

We dedicate this section to illustrate the performance of (our implementation 
of) Montes' algorithm with several polynomials chosen to force its capabilities at 
maximum in three directions: polynomials with a unique associate type of large 
order, polynomials which require a lot of refinements, and polynomials with many 
different types. 

All the tests have been done in a Linux server, with two Intel Quad Core pro- 
cessors, running at 3.0 Ghz, with 32Gb of RAM memory. 
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Example 1: Take p = 2. Consider the irreducible polynomials 

(j)^= + 4a; + 16; 
(1)2= (1)1 + 16x01 + 1024; 

03= 0i + 2"m(/.2 + 2i*^a;0i; 

04= 01 + 225x03+2350102; 

05= 01 + 2290301 + 213903 + 2^5302; 

06 = 0i + 21410305 + 227904; 

07 = 03 + 299801+21003. 

08 = 02 + 21505(^5 + 2167)06; 

09 = 02 + (((2683(xz;02 + 2l3u;)03 + 27in(u;02 + 2llxz;))0l + 

2743(a;(02 + 27t;)03 + 225(^02 + 27(«0i + 64)))04+ 

2793(^0102 + 213^)03 + 226(a;t;02 + 2l30i))05 + 2867(^03 + 21O(m;02 + 2wu))02 + 

29O5(((„0i + 64)02 + 2l2w0i)03 + 231^(02 + 270i))04 + 

296O((a;0i + 64)02 + 2i2a;i;)03 + 2986((u0i + 32a;)02 + 2i3(0i + 16«)))0708+ 

(((23364(„(^2 + 212)03 + 222(t;02 + 2llw0i))02 + 

2342O((02 + 64?iw)03 + 22O(m0i + 32a;)02 + 233^)04 + 

23469(^0102 + 2l3w;)03 + 23495((a;0i + 32u)02 + 2l2(u0i + 32x)))05 + 

2353i((y</,2 + 27(a;0i + 32w))03 + 222(0i + 16«)(02 + 2i2))02+ 

23582((t;02 + 2l6,x)03 + 225uW02)04 + 

23641((w<^2 + 213)03 + 22i(xi;02 + 2i3(0i + 16w))))06)); 

where m = x + 2,w = 0i+ 32, = 0i + 16.t. 

For each j, the corresponding polynomial (pj has a unique associate complete 
type of order j, so that in the corresponding number field Kj the ideal 2Zkj is 
the power of a uniqTic prime ideal pj. The following table contains the running 
time of the computation of a 2-integral basis in these fields with our MAGMA 
implementation of the algorithm, compared to the routines incorporated in PARI, 
MAGMA and SAGE. 

In the column "2-basis" , we give the time spent in the computation of the prime 
ideal decomposition of 2 and a 2-integral basis. The column "2-stem", shows the 
total time of the algorithm, that adds to the previous data the computation of a 
2-stem. In the columns PARI, MAGMA and SAGE we display the time spent in 
the computation of a 2-stem (by defect these programs compute always a triangular 
2-integral basis). 

Times are expressed in seconds. An "M" in a column means that the program 
exhausted the computer memory giving no result. Those executions lasting more 
than three days were interrupted. 



(j>j 


deg 0j 


ind(0j) 


2 -basis 


2 -stem 


PARI 2.3.4 


MAGMA 2. 11 


SAGE 3.2.3 


02 


4 


12 


0.00 


0.01 


0.00 


0.01 


0.01 


03 


8 


72 


0.00 


0.01 


0,004 


0.02 


0.01 


04 


16 


352 


0.00 


0.02 


0,016 


4.67 


0.05 


05 


48 


3696 


0.03 


0.6 


2.4 


42747 


4.06 


06 


96 


15408 


0.08 


0.38 


101 


> 72h 


196 


07 


288 


142416 


0.97 


16 


47157 


> 72h 


119047 


08 


576 


573696 


6.8 


M 


> 72/i 


> 72/7, 


> 72/j 


09 


1152 


23U3520 


34,5 


M 


> 72 


> 72// 


> 72/i 
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Example 2: Let f''{x) = {x^ + x + 1)^ — p^*^"*"^, with p = 1 (mod 3) a prime 
number. When we apply Montes'algorithm to factor the ideal pZk, we obtain two 
types of order zero with liftings (j)i{x) € Z[x] of degree one. For both of them the 
Newton polygon has only one side, with slope —1 and end points (2,0) and (0,2), 
and the residual polynomial is the square of a linear factor. After approximately 
2k total refinements, both types become /'^-complete. The ideal pZx splits as the 
product of two prime ideals with ramification index 2 and residual degree 1, and 
the p-index of /'^(x) is 2k. 

This is almost the illest-conditioned quartic polynomial for the algorithm, since 
the index of every type is increased a unit per refinement in general, and the total 
p-index of f''{x) is 2k. Thus, the program has to make about 2k iterations of 
the main loop. Numerical experimentation shows that even in this worst case the 
running time of the algorithm is very low. In the following table we show the 
running time of the programm for different values of k and p, compared to those 
of PARI, MAGMA and SAGE. In this case the computation of the triangular basis 
with our algorithm is almost costless, so that we have included a single column for 
the running- times. 



p 


ind(/'=) 


p-stem 


PARI 2.3.4 


MAGMA 2. 11 


SAGE 3.2.3 


7 


1000 


0.41 


2.14 


0.89 


2.4 


7 


2000 


1.14 


15.03 


3.35 


16.4 


7 


4000 


4.02 


111.7 


15.6 


121 


7 


8000 


18.9 


747 


84.6 


841 


7 


16000 


105 


5573 


486 


6374 


7 


20000 


187 


11520 


859 


12242 


13 


1000 


0.5 


3.8 


1.37 


4.4 


13 


2000 


1.5 


27.4 


5.16 


30.7 


13 


10000 


53.7 


2585 


231 


3071 


19 


10000 


65.7 


3444 


284 


4213 


31 


lOOOO 


86.5 


4741 


304 


(iOOO 


37 


10000 


93.7 


5238 


395 


6715 


43 


10000 


100.6 


5689 


422 


7370 


103 


10000 


140 


9120 


596 


11913 


1009 


1000 


0.99 


27.9 


3.65 


37 


1009 


2000 


4.49 


189 


19.6 


266.2 


1009 


4000 


24.5 


1380 


112 


2032 


10^^ + 9 


1000 


3.94 


188 


23.2 


519 


10^ + 9 


2000 


22.9 


1409 


133 


4085 


10^ + 9 


4000 


139 


10608 


763 


42790 


10«-^ + 9 


100 


1.59 


12.4 


5.61 


165 


lO^''' + 9 


200 


4.14 


88.5 


30.1 


1322 


10*^*^ + 9 


400 


14.3 


688 


167 


10802 



Example 3: Take p = 3 and let ipi{y), . . . ,i>36{y) G ^siv] be all the irreducible 
quadratic polynomials over the finite field of nine elements. Consider the order two 

types 

t'= = {x; -3, x'' + 3^ -2/3, Mv)), k = l,...,36, 
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and take representatives 4>i{x), . . . ,(j)3e{x) of these types. Define finally the poly- 
nomial 

36 

fix) :=32O"o + n0fc(x)^ 
fc=i 

The degree of f{x) is 864. Let K the number field defined by f{x). The 3-valuation 
of the index of f{x) is 757296. Our implementation of the algorithm computes this 
index and a 3-basis for the ring of integers Zx in 21 seconds. The triangulation of 
the basis, obtained after a Groebner basis computation, lasts about 3553 seconds. 
The total time for computing the 3-stem of is 3576 seconds. 



6. Some statistics 

We have tested our program with different sets of random polynomials, both 
to check the validity of conjecture 14.31 and to estimate its performance. Results 
are gathered in the following tables. In the first table, we have chosen random 
polynomials of types of orders 3,4 and 5 at p = 2, while in the second table we have 
worked with random polynomials of types of order 1,2 and 3 at different random 
primes less than 1024. 



p = 2 



Order 


Tests 


Mean Degree 


Mean Index 


Mean Time 


3 


1800 


65 


6735 


1.065 


4 


5054 


117 


25774 


3.936 


5 


300 


172 


67411 


19.605 



l<p< 1024 



Order 


Tests 


Mean Degree 


Mean Index 


Mean Time 


1 


20000 


7 


33 


0.002 


2 


10000 


25 


777 


0.151 


3 


6000 


65 


6605 


4.09 
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